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ABSTRACT 

This paper investigates the aspects of dispersion of nanoparticles in blood capillaries linked to 
nanodrug delivery systems. Nanoparticles are injected intravenously in blood capillaries and the 
resultant fluid has been identified as nanofluid. The viscosity of nanofluid is modelled using the 
nanoparticle size dependent viscosity relation. The properties of blood are studied using power law 
fluid, owing to their physically close simulations. The dispersion model framed here has been solved 
applying the method used by Sankarsubramanian and Gill for deriving exchange coefficients. The 
study conducted gives an insight into temperature, velocity and three transport coefficients of 
nanoparticles dispersed in blood with respect to various parameters like heat source parameter, 
volume fraction, power law index, size of nanoparticles, Grashof number, Darcy number and slip 
parameter for small values of wall absorption parameter under steady state conditions. MATLAB 
software has been used to plot the graphs. The outcomes reveal that physical properties of 
nanoparticles like size chiefly govern their dispersion. Convergence analysis is also stated for the 
inhomogeneous Bessel differential equation obtained while solving the problem. The developed 
mathematical model has useful applications to understand the dispersion of nanodrugs in the 
treatment of cardio vascular diseases. 

Keywords-Nanoparticle, Nanoparticle size, Nanofluid, Power-Law, Dispersion, Transport 
coefficient 


1. Introduction 

Current breakthrough in the field of nanotechnology has a profound effect on contemporary medicine 
because of their diverse merits. Enhanced bio-availability, greater diffusivity and controlled targeted 
release are some major benefits. These therapeutic nanoparticles are usually injected or inoculated 
into the patient’s body. Blood is the first medium of contiguity as the nanoparticles enter 
intravenously. The study of dispersion is an important determinant to gain knowledge of the 
functional input requirements. The nanoparticles hit the target site only after dispersing through 
blood. Hence, it is exceedingly worthwhile to realize the nanoparticles dispersing in blood. This 
mathematical model estimates the dispersion exchange coefficients which is significant to know how 
much dispersion occurs. 

The study of dispersion of any solute in blood has been an active domain of research. Implementing 
the indicator dilution method in which an amount of solute is administered in blood to comprehend 
the fluctuations in blood concentrations. Taylor [1] experimentally and analytically laid out the theory 
of dispersion of solute in steady and laminar flow through a tube. Aris [2] extended the Taylor’s 
theory of dispersion for a steady flow manoeuvring method of moments. A more generalized 
dispersion model for solute dispersion in a tube for fully developed flow was given by Gill and 
Sankarsubramanian [3]. Gill and Sankarsubramanian [4] probed into dispersion of solute in non- 
uniform fully developed laminar flow. Sankarsubramanian and Gill [5] devised a new-fangled 
approach to investigate the dispersion by resolving the boundary value problem for unsteady 
convective diffusion. Sankarsubramanian and Gill [6] extended their prementioned work on 
dispersion taking into account interphase mass transfer through the boundary of the tube wall by 
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assimilating three transport coefficients. This dispersion model has substantial relevance in diffusion 
models for blood flow. 

The dispersion theory by Taylor-Aris and dispersion model developed by Gill and 
Sankarsubramanian [7] has been executed by many researchers to interpret nanoparticle dispersion 
in various Newtonian and non-Newtonian fluids. Decuzzi et al [8] explored the diffusion of 
nanovectors in non-permeable as well as permeable capillaries employing Taylor's theory of 
dispersion. They harped upon the fact that the size of nanovector depended on the type of malignance 
of the tumour. Gentile et al [9] evaluated the nature of blood with nanoparticles in blood vessels 
modelling blood as a Casson fluid. They manifested that diffusion of nanoparticles decreased with 
the increase in the permeability of the vessel wall whereas vice versa results were obtained for 
pressure and rheological parameter. Gentile and Decuzzi [10] analysed dispersion of nanoparticles in 
blood vessels using Casson model with time dependency. 

Shaw et al [11] presented characteristics of nanoparticle dispersion in blood flowing through a micro 
vessel applying two-phase flow model of blood. The core region consisting of densified red blood 
cells was modelled as Casson fluid and the cell free region was described using Newtonian model. 
Bali et al [12] analysed the dispersion of nanoparticles in blood using Herschel-Bulkley fluid in a 
permeable capillary. Reddy et al [13] reported nanoparticle transport in blood vessels with stenosis 
using couple stress fluid model. The dispersion of nanoparticles was found to be enhanced in the 
regions of stenosis. Bali and Gupta [14] studied nanoparticle dispersion using K-L model in micro 
blood vessels with stenosis and concluded that maximum diffusion occurred at the vessel wall. 
Surabhi et al [15] examined the transport of silver nanoparticles in blood vessels applying micro- 
polar fluid model. Homotopy perturbation method was used to analyse the dynamics. Reddy and 
Srikanth [16] studied the transport of titanium dioxide nanoparticles in blood using couple stress fluid. 
Rathore and Srikanth [17] analysed the transport of nanoparticles in diseased catheterized artery. 
These outcomes emphasis on the significance of nanodrugs for treating blockages in blood vessels. 
When nanoparticles are suspended in a base fluid (in this case blood), it introduces an interaction 
force due to the kinematics of the new fluid formed, known as nanofluid. Choi [18] first proposed 
this term. Buongiaro [19] observed that the characteristics of nanofluids like dispersion was much 
enhanced and upgraded as compared to their base fluids. In this research analysis we have considered 
blood as the base fluid with nanoparticles dispersed in it, the properties of which will be described 
using the properties of nanofluids. The viscosity variation of nanoparticle dispersed in a base fluid, 
or nanofluid, was first obtained theoretically by Einstein in 1906 [20]. This model underwent 
modifications over time by various researchers. 

Non-Newtonian nanofluids are those that do not obey Newton's law of fluids. It has been studied by 
many authors to describe transport of nanoparticles in base fluids. Power-law nanofluid is an example 
of such kind where the properties of the base fluid are characterized using power law. Santra et al 
[21] studied copper-water nanofluid flowing through a horizontal channel. They used power law 
model in describing non-Newtonian nanofluids. Hayat et al [22] also studied the power law nanofluids 
inculcating the effects of convective boundary conditions too. Khan and Khan [23] described the 
MHD flow of power law nanofluid in non-linear stretching of flat surface. The properties and 
importance of power law nanofluids was studied and highlighted by many researchers [24] [25]. 
Characterizing the size of nanoparticles dispersed in blood is important as the fundamental properties 
of nanoparticles are a function of their sizes [26]. Sizes define delivery of drugs stability of the 
nanofluid, catalytic reactivity and several other physiochemical properties. Ankamwar [27] stated the 
fact that selection of nanoparticles is indispensable step prior to its usage. Size of nanoparticle has 
potential effects in their applications. Size is the representation of surface to volume ratio which 
directly influences the drug loading capacity of nanoparticles in blood circulation. Nanoparticles of 
size around 5nm have been found to be highly effective in their performance over a controlled time 
interval. Ganguly and Chakraborty [28] showed the dependence of effective viscosity model for 
nanofluid on the diameter of the dispersed nanoparticle and their volume fraction. Pasol and 
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Feullibois [29] gave the empirical relation of viscosity of nanofluid upp for its viscosity dependence 
on their diameter and volume fraction as- 


d d 
Ung = U1 + 2.5(1 — 3.4606 — + 8.6065(—52)9] 
dpy dp; 


Where yup is viscosity of base fluid in which nanoparticles are dispersed, @ is nanoparticle volume 
fraction, dnp is diameter of the spherical shaped nanoparticles dispersed and dp, is diameter of the 
blood vessel under consideration. Shaw [30] applied this viscosity model in drug delivery of 
nanoparticles in a micro vessel. 

Numerous studies have proved that nanofluids hold exceptional thermophysical properties in 
comparison to their base fluids. Hence, it is of utmost importance to inspect the parameters that effect 
thermal conductivity of nanofluids. The basic model for thermal conductivity of nanofluids was 
devised by Maxwell [31]. This model accounted for the dispersion of spherical particles in a base 
fluid for considerably small volume fraction of solutes. Bruggeman [32] gave another model to 
discern thermal conductivity of nanofluids Kk, for large volume fraction of spherical nanoparticles 


as 
kp kp kp 
(3-1) 72+ [3-) - 1+ (G6 -Dj2*[GQ - 9) - Ip +82 
Rap FN 


k 4 
f 
Where kp is thermal conductivity of base fluid in which nanoparticles are dispersed, kp is thermal 
conductivity of nanoparticles and @ is nanoparticle volume fraction. This model has higher accuracy 
as it allows for the interactivity of the spherical nanoparticles as well. 


The theoretical model for specific heat capacity of nanofluid was first given by Pak et al [33] for 
dilute solutions. An improved model for large volume fractions was given by Xuan et al [34] as 


B e! E PIP Sop + PPpcp, 


Paf — (1— $)p; + PPp 
Where Cp, is specific heat capacity of base fluid in which nanoparticles are dispersed, Cp, is specific 


C 


heat capacity of nanoparticles, pp is density of base fluid, pp is density of nanoparticles and $ is 
nanoparticle volume fraction. 


Our purpose is to investigate nanoparticle dispersion effects in blood using power law nanofluid 
model. Usually, the intravenous medications are therapeutic at small concentrations, but turn rather 
toxic at higher concentrations. Hence it is extremely important and useful to understand the role of 
dispersion of nanoparticles in blood. The power law nanofluid model has been employed to model 
the nanoparticles in blood. The equation of continuity, Navier-Stoke's equation and diffusion 
equation for temperature and concentration has been used for the mathematical formulation. The 
effects of parameters like nanoparticle size, nanoparticle volume fraction and power law index have 
been investigated on temperature, velocity and concentration of nanoparticles in blood. The 
convergence analysis for the solution of inhomogeneous Bessel's differential equation is given since 
the convergence of power series does not guarantee the existence of the solution. The results of this 
mathematical analysis will further contribute to the novel comprehension and perception of how the 
nanoparticles are transported in blood which will dispense revelation about the fabrication of 
nanoparticles in their targeted drug delivery in blood. 
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2. Mathematical Formulation and methodology 
We consider blood flow through a capillary of length L' with radius Ro as laminar, steady and 
incompressible. Cylindrical co-ordinates (r', 9’, z’) are taken into consideration for describing the 
velocity of nanofluid. In the axial direction i.e., along z'-axis, u' describes the axial velocity, while 
v’ describes the radial velocity in the capillary. The velocity along the 9’-direction is zero due to axis- 
symmetricity. 1L; describes the viscosity of nanofluid while up is viscosity of blood. 

Fig.1. Geometrical representation 


"uu 
, 


u’ 


The governing equations are given as- 


Equation of continuity- 


dp 10(rparv)) | 10pgrw! , Opgru! 
nf I T SU Burg VIE aN) y PPn un (1) 
Qt! ri Or! ri Q0! Oz! 
Navier-Stokes equation- 
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where F” in various indices stands for body forces in different co-ordinates and p,; is density of 
nanofluid. 
Diffusion equation for temperature of nanofluid- 


aT’ aT’ k ƏT! | 10T' , 0?T! H 
esu ye cnn tati a ts (5) 
Pnf Png Pnf Pnf 
where Cp, y 18 specific heat capacity of nanofluid, knp is thermal conductivity of nanofluid and ppp is 


density of the nanofluid. H is constant heat generation or absorption parameter. 
Diffusion equation for nanoparticles in blood- 

1 óc! ; 0c! ; 0c! j,O0c' co , 1 Oc! 1 0?cr , O?ci i 

Dr Atr 1 1 Wn ona T ap pe DAB 72 tm 

Dı Ot! Oz or 00! or ri Orr rr^ QO! Oz 

where D' is diffusivity and m' is rate of increase or decrease of nanoparticles. 
Nanofluids are ingenuous colloidal fluids obtained by the dispersion of 1-100nm nanoparticles in 
standard fluids. The treatment of viscosity variation in a nanofluid can be considered similar to the 


(6) 
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effect of in the viscosity of a solvent by adding particles of a solute. The relation of nanoparticle 
diameter D, with the viscosity of nanoparticles in the nanofluid unf applied is given as 


2 
Un = Mp (1 + 2.5 (1 — 3.4606 (a) + 8.6065 (25) )o (7) 
Where Do is diameter of blood capillary. uy is viscosity of base fluid. 


The thermophysical properties largely dominate the functioning of nanoparticles in the form of 
nanofluids. The chief thermodynamic properties accounted for in our mathematical model are thermal 
conductivity and specific heat capacity. The model used for describing the thermal conductivity of 
nanofluid k; is given by the relation 


. Gé-) «0-41 m DEHGA-)-1 +8 
--——— ———— 
> - (8) 
Where ky is thermal conductivity of base fluid in which nanoparticles are dispersed, kp is thermal 
conductivity of nanoparticles and @ is nanoparticle volume fraction. 


The model to determine the specific heat capacity of nanofluids Cony is described as 
(1-$)p fey + ppc 
c= Pf PD, (9) 
nf (1-$)p¢+Opp 
Where Cp, is specific heat capacity of base fluid in which nanoparticles are dispersed, Cp, is specific 


heat capacity of nanoparticles, py is density of the base fluid, p, is density of nanoparticles and @ is 
nanoparticle volume fraction. 


The nature of base fluid of a nanofluid largely governs how the nanofluid behaves. Thus, modelling 
of base fluid is highly important to predict the behavior of nanofluid. Blood is the base fluid in our 
model, the properties of which are described by a non-Newtonian fluid called the power-law fluid. 
Power-law fluids reasonably describe the blood flow behavior in tubes of diameter less than 0.2 mm 
with shear rates less than 20 s“. 

The governing equations (1) - (6) will be solved under the following assumptions:- 

1. Flow is considered two dimensional. 

. Flow is steady in capillaries. 

. Flow is axisymmetric. 

. The azimuthal component of fluid velocity is zero. 

. The cross-section area is very small in capillaries; thus, the flow is described by low Reynolds 
number. 

6. Power-law nanofluid model describes the properties of fluid in capillary. 

7. A dilute solution is assumed to solve the dispersion model. 

8. Free convection effects are ignored. 


Un ULP 


The modified temperature diffusion equation and Navier-Stokes equations along with their respective 
boundary conditions are described henceforth. The dispersion of nanoparticles in blood capillary is 
described later in the section under dispersion model. 


The modified temperature diffusion equation followed by its boundary conditions is given as:- 
, 0T! , OT! knf 9?T! 1ôT' on (10) 


ii arr aoe az Pnf Pap Gre riarr | Oz? 
The temperature is prescribed Tọ all over surface and the temperature gradient disappears throughout 
the capillary axis. 
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OT! " 
5;,70 atr'=0 (11) 
T =T, atr'= Ro (12) 


The modified Navier-Stokes equation followed by its boundary conditions is given as:- 
ð d ð II 1 
Pug (1) * ona T — To) = 0 (13) 


Oz! r'ar' 
The shear stress is constant along the axis of the capillary because there is a constant source causing 
a deformation in the nanofluid, the boundary condition for is expressed as:- 

T’ is finite at r' = 0 (14) 


The shear stress relation for power law fluids is given as:- 
Qu! 


Uu n (15) 
Where n is the power law index. 


Using Darcy's law at the boundary of the capillary, we get 


u'-ugatr'-Hg (16) 

our o! 1 1 ri 

ar 7 pg 8 Up) at r' = Ro (17) 

where 

gm Dco is velocity at the permeable boundary (18) 
Ung Oz! 


Where u'g is slip velocity, ø’ is slip parameter, Da is Darcy number. 


The non-dimensional scheme is stated as:- 


pet nat ga" Ga Oe x LL HN e = taf plu uc". — Ror! 
Ro Ro T Uavglf Ty kf Uf 0 Uavg Uflavg 
ol 
^ (19) 
The non-dimensional equations are stated as:- 
k k k 
(36-15, «[3a-6)-1]* (36-1), [30-9)-1]? +82 
3290 | 100 Sr f LS 20 
mh HÉROES 2 
EXE " 2 A e 2 
Re = 23 (rur (1 25(1 3.4606 (È) + 8.6065 (2) )e) G9 )+(a p) + 
seer) Gro (21) 
DfYf 


The non-dimensional boundary conditions are stated as:- 


ae _ 
—=0 atr=0 (22) 
0—1atr-1 (23) 
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u-—upgatr-1 (24) 
ðu o 
ar^ qae Up) at r=1 (25) 
where 

— _ Dadp = 
Up = "o atr=1 (26) 


2.1 Solution for temperature 0 


The equation (20) is solved analytically using the boundary conditions (22) and (23), and the value 
of 0 is obtained as:- 


h(1-r2)(36-1) 22 «[3(1-$)-1]* |((36-1):2[(3(1-9)-1]2*825) 
kf kf kf 


0 = +1 (27) 


4 
2.2 Solution for velocity u 


The value of 0 is substituted in equation (21) and is solved analytically using the boundary conditions 
(24) to (26) to obtain the value of u and ug as:- 


$ppY k k k 
h(-6), 5 are(6-1,7 «0-91 jen Dpp GG )-1]}? +87") ; [c = 
er 5 
( 164g (1*2. s(1- 3460e( P P) +8. 6065(7 z) Jo) ) ubt mre) 
à (a -9)4 = Gro + 
PfYf 


m(a- store as- Dee a-4)-1 jen D «(Ga-) 7 


1 


va- ost pe aso Dg a-$)-1* jen DEGU) zo 


dz 16y(142.5(1- 3.4606 (P )+8.6065 (7 zb) e) nt 


A) cepe? 


Ha pa 15? «[0-9)-1]* Jes 15 -[G-49) Mee 


1 
xa- ot p aso Dg«a-$)-1 [ce Di «Ga-) wd 
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Re 3 Lous (142. 5(1- 3.4606 (P )+8.6065 (7 2)’)6) (suci) 
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Ha peta a oe Dire BC- à)-1]* jen 1) ea )-1)}2+8 =) 


E) Fa 
16up(1+2. s(- 3.4606(2 P)+8. eoes (7! i) Jo) ganz B 
(28) 

Applying boundary condition (25) and (26) to find the value of ug 

Atr=1in 

ðu 

ðr — 
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(a- D — 


2 2 0z Jal B uy OZ 
We get up as 
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2.3 Dispersion model 


The modified d equation (5) using the assumptions and the boundary conditions are given as 
dcr Och ð AG E Q?c! 

Qt! LE Oz DC. ür! ür! + 022) (32) 
Where t'denotes the time in dimensional form. 


The nanoparticles are distributed uniformly in the capillary at z' -— 0 therefore the initial 
concentration distribution is described as 


c'(0,z', r^) = w'(z)&(r) (33) 
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Where w'(z’) = AS 6(z’) is the Dirac Delta function (34) 
1 o <r<d' 


Where d' is the initial distribution. 


Initially, at the beginning of the dispersion process 


c'(t', z', r^) = finiteatr' = 0 (36) 
Diffusion of nanoparticles at the capillary wall can be expressed as 
-Dp 2 27) — ky'e'(t' z', 7) atr! = Ro (37) 


Where kọ’ is the permeability at the capillary wall 


The expression for finite quantity of nanoparticles in the blood stream at any instant of time is 
à 
c'(t'z',r) = sete ar) = Qatz' =% (38) 


The non-dimensional scheme is stated as:- 
Dmt! Rou cI kioR ul T! zl 
t= Pe = = e= p= u= „r =—,z =^ (39) 
R2 Dm Co Dm Uavg Ro Ro 
Where is Pe Peclet number, c, is reference concentration and f is non-dimensional wall absorption 


parameter. 


The non-dimensional equation is written as 


Oc Oc 10 Oc 1 0?c 
a es bac 9) 
The non-dimensional initial and boundary conditions are given as 
c = w(z)é(r) where w(z) = TE (41) 
|» (10«rxd 

And £9) = {9 q <r<1 
c = finite atr = 0 (42) 
26 2 -Bcatr-1 43) 
rm Bcatr = ( 

= = 0 atz = co 44) 
c—3,- 0atz- ( 


2.3.1 Solution of the dispersion model 


The solution of convective diffusion equation in non-dimensional form is given using the method 
B by Sankarsubramanian and Gill [6] as 


0c. 
i-o Wi 5, T 


. ; : : 1 : 
where the MT mean concentration Cn is described as Cm = 2 A cr dr and yj(t,r),i = 


(45) 


0,1,2, ... is a function of time and radial distance. (46) 


The governing equation (40) is multiplied by 2r and integrated with respect to r from 0 to 1 
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OCm _ 1 02cm ðc 


ðt Pe? az? ðrlr=1 


a pl 
— 2z h ucr dr (47) 


The dispersion model for Cm with time dependent exchange coefficient using equation (45) in 
equation wae 


Si 
D (48) 
T 
K= +2% -2f pimurdr,i= 0,1,2 p =0 (49) 
i Pez ni 0 i-1 , — yp hy Gy van wenn , -L— 
i i=j 
Here 6;; " LEJ (50) 
The resultant dispersion model for mean concentration is given as 
3cm 3cm 32cm 
a = Kot Ky + EST (51) 
oi 
K,2329 429^| | .2[ yi ilurdr,i 20,12 (52) 
a ar lr=1 0 i-1 , yd) 


The exchange coefficient Ky is due to non-zero nanoparticle flux at the capillary wall, K, exists due 
to the convection coefficient because of the velocity of nanoparticle and K, is the dispersion 
coefficient due to molecular diffusion and velocity of nanofluid. The terms K; for i = 3,4,5, ..... have 
been neglected as the contribution of higher order terms is insignificant. 


The d d concentration is thus expressed as 


alc 
i=0 Wis T 


(53) 


For solving equation (51) and (52), we peo to find the boundary conditions. Using equation (45) in 


equation (40). Equating the coefficients of 2 — for | = 0,1,2 and obtaining 


api 10f( 0 
Minie 2 (r ot) EL PETI i p where l = 0,1,2 (54) 
and po = Y2 = 0 (55) 
The initial and boundary conditions are obtained from (41) to (44), described as 

m = 20 f &r dr att =0 (56) 
Po =—— att =0 (57) 

974 fo ér dr 
y, = 0,1=1,2 att =0 (58) 
“= 0,1 = 0,1,2 atr = 0 (59) 
St = fy, 1 =01,2atr =1 (60) 

3cm 

Cm ==, — 0atz — oo (61) 


Using condition (55) in (52), we have an additional condition as 
f; Wir dr =} ôw for l = 0,1,2 (62) 
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Where ô; is defined by (50) 
2.3.1.1 Solution for Wo and Ko 


The value of Wo and Ko are independent of velocity of nanofluid, thus they can be obtained directly 
using the boundary conditions (55) 


avo 19 (,,2v0) _ 
(Qt  rór (r or ) Koo i (63) 
With an additional condition on Wo as f. o Vor dr — 1/2 (64) 


The solution on non-homogeneous boundary value problem (63) using Bessel equation satisfying the 
conditions (56)-(62) and (64) is given as 


Let polt, r) = eC fo Koman oct, r) (65) 


Using this transformation (65), we get an equation in terms of o(t,r) that has to satisfy 


3e 10 „2e 
at ^ rÓr T or (66) 
g(0,r) = vy (0,r) = $2 — (67) 
, Ust 2 fo r£(r)dr 
8 
5. (t,1) = —fo(t, 1) (68) 
o(t, 0) = finite (69) 
Let the solution be of the form 
elt, r) = Xm-o Ame ^nt (70) 
Using the well-known [35] identity 
do Or. —— or 
GerGe t= ape (71) 
We convert (66) into (72) as 
10r? _ ð , 8g 
Integrating it with respect to o yields 
Ec do = ew (73) 


By substituting (70) into (73), we obtain the solution for o(t, r) by equating the coefficients on each 
side. Thus, the solution of (66) to (69) is given as 


p(t,r)e X5sA E tnt (74) 


Now, using (65), we get 

E EL o AmJo(Amr)e ^nt 
Vo = 2 Yo -o (52.4 n) "nt 
and Àm are the roots of Bessel equation 
Amim) = Bigs) ot = 0,1,2, n. (76) 


_ Aha fo TJo(Amr)dr B 
Also Am = OR 82) Am) fd rt ar’ m= 0, 1,2, does (77) 


where Jọ and J, are Bessel functions (75) 
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Kg is obtained from the initial condition as 
Apo YS o AmÁn fi (Ag)e "ht 
Ko = 2— = Um LÁ 78 
NEL. MEN NE uo 
For steady state, t > oo , thus, equations (49) and (65) give the undermentioned asymptotic 
representation for the function Yo and Ko 
|. Ào EM 

lim Uo = Iag (Aor) and lim Ky = —4$ (79) 
where A, is the smallest root of the Bessel's equation (76) 


2.3.1.2 Solution for K; and y; , i=1,2 


Under steady flow conditions, the value of K; , i=1,2 is derived applying the value of nanofluid 
velocity obtained using power-law nanofluid model. The function yj, i=1,2 for steady state condition 
is found as 


10 ð 1 

12 (r) + adus = uia — zz Yia + Dla Kipi- + Kipo where 151,2 and p_,0 (80) 

For steady state K; , i=1,2 reduce to 

K=% +2! —2f rw dr , 1-212 (81) 


^ pe? ðr 


r=1 


The boundary conditions on y;, l= 1,2 


y; = finite or = Oatr = 0,1212 (82) 
TL By atr=1, 1=1,2 (83) 
f, wr dr =0, 1212 (84) 


The equation (80) describes the Strum-Liouville boundary value problem, the solution of which is 
found using the property of orthogonality of characteristics functions. The characteristic function of 
the corresponding homogeneous boundary value problem will be orthogonal to the right-hand side of 
the equation (80) with the weighting function r. So, we find the solution of homogeneous differential 
equation corresponding to (80) i.e. 


12 (Ot) 492 = 

rr (r or ) 7 oyi = (85) 
The value of ij; from here can be directly calculated using the Bessel’s function of zeroth order as 
V, = Jo(Aor). 

Let this be denoted by characteristic function $,, = Jo(Aor) (86) 


Now, multiplying both the sides of the equation (80) by røn, the left-hand side of the equation 
vanishes and using the property of orthogonal of functions, we obtain the value of K,. The expression 
for the exchange coefficients Kj, l = 1,2 is obtained in terms of the function y,,/ = 1,2 

|o Toor) (ubi-1-poz1-2+ Dia Kiar 


K, = yd 12 87 
i [is VorJo(Aor)dr (87) 
The value of K, using (84), (87) in (81) 
_ fo rJoAorupodr 2 i. 25 
Ki = rnana REDA Jo urloAor)dr (38) 
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The value of y; from (70) 
0 ð 
19 (291). yp, = upo + Kipo (89) 


ror or 


The boundary condition for Yı 


y, = finite atr = 0 (90) 
Oi = 
p = BY atr=1 (91) 
f; Yur dr = 0 (92) 
Using the value of K, in (87), we get the solution of y, satisfying (90) to (92) as 
V; = Xm-oBmlo Amr) (93) 
where the expansion coefficient By using (92) and (93) is obtained in terms of B; (j=1,2, .....) as 
90 7 Jao) nim Am 


Using (94) in (93) 


co Ao Js Am) 
V1 = Xm-o Bm Uo Amr) — D p Jo(4or)) (95) 
24A 1 
Where Bm = xar A (u + Ky) WorJo Amr)dr (96) 


The dispersion coefficient K, using (88), (95) and (96) is 
1 hlo) 1 
Ko = xz Gp Faw Jo © + KDY Aor)dr (97) 


The exchange coefficient Ky occurring because of the absence of nanoparticle at the capillary wall. 
The convection coefficient K, and dispersion coefficient K, are obtained for fj # 0, hence they are 
responsible for the phenomenon of dispersion of nanoparticles in the capillary and the values for 
which are obtained as 


Kı = -2 fru dr (98) 


T 
n 
2 hB 18n?+n*4+18n+24n?+6n3 16 
n- | emne e(l gy 4 
: (om) (9-n?)(3+4n+n?) hB ( 
1 


PPpYp Op hB ( 9n?—7n^*-18n? ) 
——JGr8 + Re— || — —————— ——-J]-—A 
DfYf ) 3 (ae) (9-n?)(9-7n? -2n?) 


(99) 


Now, K, = =; -2 f, ruypdr (100) 
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hB n? A? n? 


= —- 2| —— o BA erc 
Pe? (oe 3.4606(2 P)+8.6065(7> za) Eia i^ Sn ron) 


Aim? u 2n Aen 7 Ain )'z _ 57 
64(3-n)(3-7n) (3—n)(3+n) Tem 32(3-n)(345n) ü $) Pe rr Ioram 


n 
ap hB 5 n? Af? 
mD"—L g(t ——— 
3 e 34606? )+8.6065(2 any) 1=1°'\(3-n)(3+3n)  4(3-n)(343n) 
Ain? 2n3 an? _ Ain? ) ( A 
64(3—n)(345n) (GU RUBEPU peu E 32(3-n)(347n)(1-n) +A XB (1 
1 
n 
FAMAE ——— E ME 
144 FES EY ROB ieu (12.5 (1- 3.4606 (2 2) +8.6065(“2) jo) (3-n)G-43n)  4G-n)G+5n) 
Aon? _ 2n Aan _ aon ) | 16 = PPpYp 
64(3—n)(347n) (Sene RDN 32(3—n)(3+5n) hB (1-9) + DfYf quitan 


1 
n 


ge? nB "ys, p, | 2200 B ( 
dz 16y(142.5(1- 3460e( P P)+8.6065(72 20 n NACE 16 (12.5 (1- 3.4606 (2 7) +8.6065(5 2)’)6) (3-1 


Aan? Aon? _ 2n3 Ažn3 _ Adr? ) 
4(3-n)(345n)  64(3-n)(347n)  (3-n)(34n)  2(3-n)(343n)  32(3-n)(345mn) 
o AoD (4 AL, A 
A Xia Bi ES (1 a 2) (101) 
Where e e (1—4$)-4 
Ha peta 1) PUO à)-1]* jen 1) Ga $) ea 
a 
Ha pote 1) epa $)-1]* jen 1) ace $)-1)?48 2 
HAVAA A S (j 
PfYf dz 16y(142.5(1- 3.4606(7 7) +8.6065(5 22) Je) (9-n?)(1-n) 


1 


Ha pote 1) PUO $)-1]* ce DEA $) zo 


16up(1+2. s(1- 3.4606 (P P)+8. 6065 (2 20 


(102) 
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B= (a — $) eere (ao - 1)2 430159) —1]4 
PfYf kf 
(86 - D? - [GQ - 9) - 1 + 2 (103) 
f f 


2.3.1.3 Solution for mean concentration 
The equation (51) is solved using the boundary conditions (56) and (61) to get 


Cm = aan — zy (104) 
Where A(t) = f; Ko(n)dn (105) 
z, (tz) = z + f, Ky(n)dn (106) 
T(t) = fr K (dn (107) 


We have considered large time interval to analyse the dispersion model. So, the approximation of 
equations (95) to (97) for large time can be represented as 


A(t) - Kot (108) 
Zy(t,z)-z + Kat (109) 
T(t)~Kzt (110) 


2.4 Convergence of solution 

Convergence plays an important role in determining the accuracy of the solution. In this problem we 
have calculated the value of y(t, r), i = 0,1,2, ... which is a function of time and radial distance. The 
value of i, is in the form of Bessel function. The convergence of the Bessel function exits in the 
form of power series solution, but, it does not guarantee the existence of solution to inhomogeneous 
Bessel differential equation [36] as in equation number (80). Using equation number (93) 


W1 = Xm-o BmJo Amr) 


Lemma 1:- If y, /r has bounded variation in the interval 0 € r < 1 then fh ViJfo(Amr) = k e 
Proof:- Using the basic properties of Bessel functions, we have 
_ 72 1/2 Rm gr) 
Jo(r) = C) [cos G z) tk } (111) 
Let p = oe , then 
f Wo Amndr = (S) for cosAmr — Ddr + k £z (112) 
0 170 Vim = TÀm 0 m 4 43/2 


9 will be monotonically increasing if 9 = 94(r) — V(r) in which 94(r) and 9;(r) are also 
monotonically increasing. Hence, 
! &(r)r? cos(Amr — Ddr = 8(1 + 0) fË rž cos (Amr — 2) dr + 8(1 — 0) ft r2cos(A 
f, 8r)r2 cos( mr — dr = (1+ ) f; 72 cos ( mt — =] r+9(1- ) J; ricos (Amr — 
Z) dr = x £&m (113) 


4 


Lemma 2:- In the interval 0 € r < 1, let the function i, be absolutely continuous and let yw,’ has 
bounded variation. Then, the general coefficient B,, of the series in (93) may be written as:- 

1/25. 41m | 
Qu) gD su PU) where nc {° when i # 0 


ANA ae 1 wheni = 0 
Proof:- Using the property of Bessel function we have 
1 
fo HR Amr)ar = = YF Am) +J? Am) (114) 
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Using (111) 
1 i (Am) 
fy ndr - + KS (115) 
By (115), Bm in (93) assumes the form 
1 
[Amt + KP1(Am)] Jy rior) dr (116) 
On integrating using the recurrence formula 
d 
a; AQ] = rJfa(r) (117) 
The integral (116) becomes 
1 1 pl 
So rios dr = = f; dl Amr Ii Am?) (118) 
1J1(Am 1 pl_,, 
Am) — t fo rV Jis dr (119) 
If we assume i’, = 0 and monotonically increasing, thus, = [y (r) — y4(0)] = 0 is monotonically 
increasing. It follows that in (119) the integral L f E ry JA (A r)dr has the form x x 
Now the term 91829 ° 
Am 2 mz qc k$-ky(m).m 1€ y(m) <2 (120) 
nom, 
xtn—-z+- i=0 
Where q = CM x is an integer (121) 
mt i0 
From (120) 
. : ( 
sin(A& — a = sin (mr +q- z) + y 209 (122) 
ti(-1)!"-,999 i=0 
=) gm i (123) 
k —— i#0 
m 
And (111) gives 
I Qs) = CD" CE? + eon (124) 


Theorem:- Let iJ, be a function as described in lemma 2 and in addition let the condition ó;y4 (1) = 


0 be satisfied then y4 — S,,(r) = K on in 0 € r € 1 where S,,(r) is the sum of first m regular 


terms of the series A. 

Proof:- It has been shown in Lemma 1 that the series will converge to the value of the function in any 
subinterval of 0 € r € 1. If y, (r) is continuous in this subinterval and 6;y, (1) = 0, therefore under 
the conditions of the theorem the series will converge to i4 in 0 <r < 1. Further, from lemma 2, the 


D Jo(Amr) where Jo (Amr) is uniformly bounded 


general terms of the series can be expressed as k 
Am 
gm) 


mi? 


and the remainder after m terms is written as k Thus, 4 is convergent in 0 <r < 1. 


In the above method we described a condition on the remainder after m terms, which is a stronger 
condition for the convergence of the inhomogeneous solution of the Bessel function. Thus, the 
solution Wy = 35e Bao (Agr) obtained for the equation converges. 


3. Results and discussions 

The present mathematical investigation analyses the dispersion of nanoparticles in capillaries using 
power law model, over steady state conditions. The temperature of the nanofluid has been analysed 
for different values of heat source parameter and volume fraction of nanoparticles. The effects of heat 
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source parameter, volume fraction of nanoparticles, power law index, radius of nanoparticle, Grashof 
number, Darcy number and slip parameter has been observed on velocity of nanofluid. 

We have employed the method adopted by Sankarsubramanian and Gill [6] and obtained the effects 
of parameters like heat source parameter, power law index, radius of nanoparticle and Grashof 
number on the exchange coefficient, convection coefficient and dispersion coefficient during 
nanoparticle dispersion under steady state conditions with small values of wall absorption parameter. 
The small values of wall absorption parameter signify small or no reaction rates at the healthy vessel 
wall. This would justify the fact that as the nanoparticles move towards the wall their concentration 
at the wall will not become zero for healthy cells at the wall and hence, they can disperse into the 
diseased cells which will have a higher value of absorption parameter than the healthy cells. Thus, 
the nanoparticle can be used for diseased cell targeting. 

Fig 2 displays the effects of temperature 0 of nanofluid against radial direction r for different values 
of heat source parameter h . The temperature drops with the increase in radial co-ordinates. This is 
because nanoparticles convect radially from higher temperature gradient to lower temperature 
gradient. It is seen from Fig 2 that higher the value of heat source parameter h, greater the temperature 
of the nanofluid. This is because the increase in the value of heat source parameter signifies greater 
heat production in the nanofluid [37] due to nanoparticle Brownian motion. Thus, the temperature of 
nanofluid rises with the increase in heat source parameter. 

Fig 3 depicts the effects of temperature 0 of nanofluid against radial direction r for different values 
of volume fraction $ of nanoparticles. The trend shows that higher the value of volume fraction, 
greater the temperature of the nanofluid. The volume fraction of nanoparticles represents the number 
of nanoparticles in the nanofluid. Thus, higher the value of volume fraction, greater the temperature 
owing to increase in the number of nanoparticles, hence the rise in temperature. Similar results were 
reported by Tan et al [38]. 

Fig 4 shows the graph of velocity u of nanofluid against radial direction r for different values of heat 
source parameter h . The rise in the value of heat source parameter, raises the temperature of the 
nanofluid which in turn causes a rise in the velocity of the nanofluid due to higher heat generation 
[37]. Thus, the velocity rises with the increase in heat source parameter. 

Fig 5 shows the graph of velocity u of nanofluid against radial direction r for different values of 
volume fraction $ of nanoparticles. The trend shows decrease in velocity with increase in the volume 
fraction. 

Higher value of volume fraction signifies higher number of nanoparticles in the nanofluid. The greater 
number of nanoparticles experience greater collision amongst themselves, thus decreasing the value 
of velocity of nanofluid [37]. 

Fig 6 shows the graph of velocity u of nanofluid against radial direction r for different values of power 
law index n. The trend depicts that higher value of power law index causes a decrease in the velocity 
of nanofluid. Power law index represents shear thinning fluids for values lesser than unity. This 
behaviour is observed because of the nanoparticles dispersed in the blood. Nanoparticles tend to form 
loose aggregates which is broken rapidly by increasing shear rates. This causes a reduction in their 
resistance to flow. Thus, the velocity of nanofluid decreases with the increase in the power law index 
[24]. 

Fig 7 displays the graph of velocity u of nanofluid against radial direction r for different values of 
nanoparticle radius rp. The nature of the graph shows that decrease in the value of nanoparticle radius 
increases the velocity of the nanofluid. The viscosity dependence relation given by Pasol and 
Feuillebois [29] directly relates the viscosity with the radius of nanoparticle. And viscosity is 
inversely related to velocity for shear thinning fluids as the case considered here. Nanoparticle radius 
governs the size of nanoparticle. Higher the radius, greater the size of the nanoparticle. The greater 
size of nanoparticle causes a reduction in velocity, thus the result. 

Fig 8 exhibits the graph of velocity u of nanofluid against radial direction r for different values of 
Grashof number Gr. The trend shows that the velocity increases with the increase in the value of 
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Grashof number. Grashof number is the ratio of buoyant forces to viscous forces. Buoyant forces are 
caused in the nanofluid due to the temperature gradient, as there will be no flow in the absence of it. 
Thus, higher the value of Grashof number, higher the temperature, that in turn causes a higher velocity 
[37]. 

Fig 9 displays the graph of velocity u of nanofluid against radial direction r for different values of 
Darcy number Da. The graph shows that the velocity increases with the increase in the value of Darcy 
number. Darcy number represents the permeability of the medium across the cross-sectional area of 
fluid flow. As the value of Dracy number increases, the fluid flow becomes stronger due to increase 
in the permeability of the medium [37]. Thus, increase in Darcy number causes an increase in the 
velocity of the nanofluid. 

Fig 10 shows the graph of velocity u of nanofluid against radial direction r for different values of slip 
parameter o. It is seen that an increase in the value of slip parameter causes a decrease in the velocity 
of the nanofluid. Slip parameter represents the fluid slip at the boundary of the capillary wall. The 
increase in its value causes a decrease in the skin friction and heat transfer, thus deceasing the velocity 
of the nanofluid [39]. 

Fig 11 exhibits the graph of exchange coefficient —K, against wall absorption parameter f. The 
graphs shows that exchange coefficient Kg increases with the increase in wall absorption parameter 
B. Ko signifies the reaction rate constant depending on time. As it is a steady state, Ky depends only 
on Wo. Wo describes the dispersion of nanoparticles in the capillary for first order homogeneous 
reactions. Wall absorption parameter p represents the nanoparticle concentration that builds up at the 
wall of the capillary as they are dispersed. Thus, the reaction rate constant depends directly on the 
rate at which nanoparticles are absorbed at the capillary wall [6] [40]. Hence Ky increases with the 
increase in f. 

Fig 12 exhibits the graph of convection coefficient —K, against wall absorption parameter p for 
different values of heat source parameter h. Convection coefficient enhances with the increase in wall 
absorption parameter because nanoparticles begin convections rapidly towards the capillary wall due 
to there increased interactions at the wall of capillary [40]. The graph shows that an increase in the 
value of heat source parameter h causes an increase in —K,. Since heat source parameter relies on the 
heat production due to nanoparticle collisions, an increase in its value increases the dispersion of 
nanoparticles in the fluid. Thus, the value of convection coefficient increases with the increase in the 
value of heat source parameter. 

Fig 13 shows the graph of convection coefficient —K, against wall absorption parameter p for 
different values of power law index n. It is seen that an increase in the value of power law index 
causes a decrease in convection coefficient. Power law index describes the rheology of the nanofluid. 
Nanoparticles behave like suspensions, thus nanofluid shows shear thinning properties under low 
shear rates. Since it is a shear thinning fluid, as the value of power law index decreases, nanoparticles 
gradually align themselves in the direction of increasing shear and thus impedance to flow is reduced. 
As a result, nanoparticles start convections faster towards the capillary wall [40]. Thus, the convection 
coefficient —K, increases as the value of power law index n decreases. 

Fig 14 displays the trend of convection coefficient —K, against wall absorption parameter fj for 
different values of nanoparticle radius rj. From the graph it is seen that an increase in the value of 
nanoparticle radius causes a decline in the convection coefficient. The radius of nanoparticle chiefly 
defines the size of nanoparticle dispersed in the nanofluid. Greater the size, more will be impedance 
to flow, thus slower will be the convection. Hence the result. This implies that that the physical 
properties of nanoparticles hold considerable importance in their dispersion [41]. 

Fig 15 shows the graph of convection coefficient —K, against wall absorption parameter f for 
different values of Grashof number Gr. It is seen that an increase in the value of Grashof number 
causes an increase in convection coefficient. Grashof number arises due to convection. Thus, higher 
the value of Grashof number, greater the convections. 
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Fig 16 exhibits the graph of dispersion coefficient —Kz against wall absorption parameter fj for 
different values of heat source parameter h. Dispersion coefficient decreases with the increase in wall 
absorption parameter. The velocity gradient is smaller in the central region of the capillary than near 
the capillary wall [6]. The nanoparticles tend to disperse more towards the capillary wall with the 
increasing wall absorption parameter, increasing the transverse movements towards the wall, and thus 
decreasing the axial dispersion coefficient [6]. The graph also shows that an increase in the value of 
heat source parameter causes an increase in dispersion coefficient. This is because as the value of 
heat source parameter increases the dispersion is enhanced in the capillary. 


Fig 17 exhibits the graph of dispersion coefficient —K; against wall absorption parameter f for 
different values of power law index n. It is seen that as the power law index increases, dispersion 
decreases. Power law index defines shear thinning behaviour for values less than unity. Nanoparticles 
suspended in nanofluids define the properties of a shear thinning fluid. Therefore, as the value of 
power law index decreases, the nanofluid becomes less viscous and hence the dispersion increases. 


Fig 18 displays the trend of dispersion coefficient —K; against wall absorption parameter fj for 
different values of nanoparticle radius rj. From the graph it is seen that an increase in the value of 
nanoparticle radius causes a decline in the dispersion coefficient. The nanoparticle size gives the 
measure of nanoparticle coverage which quantifies and qualifies the dispersion. Higher dispersion 
occurs with smaller nanoparticles because they cause lesser impedance [42]. Thus, greater the size, 
lesser the dispersion of nanoparticles. 


Fig 19 shows the graph of dispersion coefficient —K; against wall absorption parameter f for 
different values of Grashof number Gr. It is seen that an increase in the value of Grashof number 
causes an increase in dispersion coefficient. This is because an increase in Grashof number causes a 
greater movement of nanoparticles towards the capillary wall [43]. The buoyancy forces enhance 
dispersion, thus the result. 
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Fig. 2. Variation of temperature 0 against radius r for 
different values of heat source parameter h 
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Fig. 3. Variation of temperature 0 against radius r for 
different values of volume fraction @ 
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Fig. 4. Variation of velocity u against radius r for different 
values of heat soure parameter h 


1000.0022 
900.00224 |— — ~= 
b 
800.00224 | RS 
zx g 
È = = = volume 
S 
oo fraction=0.01 
600.00224 
3 500.00224 -——]——[—o—— eee M volume 
fraction=0.02 
400.00224 
300.00224 volume 
200.00224 fraction=0.03 
100.00224 
0.0022445 


0 01 02 03 04 05 06 07 08 09 1 


Fig. 5. Variation of velocity u against radius r for different 
values of volume fraction @ 
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Fig. 6. Variation of velocity u against radius r for different values of 
powerlaw index n 
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Fig. 8. Variation of velocity u against radius r for different values 
of Grashof number Gr 
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Fig. 12. Variation of convection coefficient -K1 against wall absorption 
parameter p for different values of heat source parameter h 
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Fig .13. Variation of convection coefficient -K1 against wall absorption 


parameter f for different values of power law index n 
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Fig .11. Variation of exchange coefficient -K0 against small values of 
wall absorption parameter f under steady state condition 
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Fig 15 Variation of convection coefficient -K1 against wall 
absorption parameter p for different values of Grashof number Gr 
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Fig. 16. Variation of dispersion coefficient -K2 against wall 
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Fig. 17. Variation of dispersion coefficient -K2 against wall 
absorption parameter p for different values of power law index n 
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4. Conclusion 

Nanoparticles have become an exceedingly powerful tool for drug delivery in the treatment pf 
cardiovascular diseases. The evaluation of standard processes in the field is a major field of research. 
This paper concerns with the nanoparticle dispersion in capillaries using power law model. The effect 
of heat source parameter, volume fraction, power law index, size of nanoparticles, Grashof number, 
Darcy number and slip parameter for small values of wall absorption parameter under steady state 
conditions have been examined on temperature, velocity, exchange coefficient, convection 
coefficient and dispersion coefficient of nanoparticles dispersed in blood capillary. The formula 
predictions agree well with the experimental designs [44]. The above nanoparticle aspects are 
important criteria for nanoparticle dispersion in physiological situations. The major findings of the 
study can be outlined as:- 

e The temperature of nanofluid rises with rise in heat source parameter and volume fraction of 
nanoparticles. 

e The velocity of nanofluid enhances with enhance in heat source parameter, Grashof number and 
Darcy number. 

e The velocity of nanofluid declines with an increase in volume fraction, power law index, 
nanoparticle radius and slip parameter. 

The exchange coefficient boosts with the increase in wall absorption parameter. 

The convection coefficient also grows with rise in wall absorption parameter. 

The convection coefficient elevates with elevation in heat source parameter and Grashof number. 
The convection coefficient decreases with increase in power law index and nanoparticle radius. 
The dispersion coefficient reduces with the increase in wall absorption parameter. 

The dispersion coefficient improves with elevation in heat source parameter and Grashof number. 
The dispersion coefficient subsides with increase in power law index and nanoparticle radius. 
This model provides an efficient insight into the physical properties of nanoparticles, most 
significantly their size, which is important for their dispersion in blood. It has useful application in 
the treatment of targeted delivery of drugs where nanoparticles can themselves serve as drugs or a 
suitable drug can be attached to it. In future, the dispersion model discussed here can be extended to 
study nanoparticle clustering and also their distribution in full vasculature. 


References 

[1]G Taylor. “Dispersion Of Soluble Matter In Solvent Flowing Slowly Through A Tube". 
Proceedings of The Royal Society London A 219, pp 186-203, 1953 DOI:10.1098/rspa.1953.0139 
[2]R Aris. “On The Dispersion Of A Solute In A Fluid Flowing Through A Tube". Proceedings of 
The Royal Society London A 235, pp 67-77, 1956 DOI:10.1098/rspa. 1956.0065 

[3]W N Gill and R Sankarsubramanian. “Exact Analysis Of Unsteady Convective Diffusion." 
Proceedings of The Royal Society London A 316, pp 341-350, 1970 DOI:10.1098/rspa.1971.0083 
[4]R Sankarsubramanian and W N Gill. “Dispersion Of Non-Uniform Slug In Time-Dependent 
Flow" Proceedings of The Royal Society London A 322, pp 101-117, 1971 
DOL:10.1098/rspa.1971.0057 

[5]R Sankarsubramanian and W N Gill. “Dispersion From A Prescribed Concentration In Time 
Variable Flow." Proceedings of The Royal Society London 329, pp 479-492, 1972 
DOL:10.1098/rspa.1972.0125 

[6]R Sankarsubramanian and W N Gill. *Unsteady Convective Diffusion With Interphase Mass 
Transfer." Proceedings of The Royal Society London A 833, pp 115-132, 1973 
DOI:10.1098/rspa.1973.0051 

[7]R Sankarsubramanian and W N Gill. “Corrections to ‘Unsteady convective diffusion with 
interphase mass transfer’ ". Proceedings of The Royal Society London A 341, pp 407-408, 1974 
DOI:10.1098/rspa.1974.0195 


(02023, IJETMS | Impact Factor Value: 5.672 | Page 454 


International Journal of Engineering Technology and Management Sciences 
Website: ijetms.in Issue: 2 Volume No.7 March - April — 2023 
DOI:10.46647/ijetms.2023.v07102.052 ISSN: 2581-4621 


[8]P Decuzzi, F Causa, M Ferrari and A P Netti. “The Effective Dispersion Of Nanovectors Within 
The Tumour Microvasculature." Annals of Biomedical Engineering 32, pp 798-802, 2006 

[9]F Gentile, M Ferrari and P Decuzzi. “The Transport Of Nanoparticles In Blood Vessels: The Effect 
Of Vessel Permeability And Blood Rheology." Annals of Biomedical Engineering 36, pp 254-261, 
2008 

[10] F Gentile and P Decuzzi. “Time Dependent Dispersion Of Nanoparticles In Blood Vessels." 
Journal of Biomedical Science and Engineering 3, pp 517-524, 2010 DOI:10.4236/jbise.2010.35072 
[11] S Shaw, S Ganguly, P Sibanda and S Chakraborty. “Dispersion Characteristics Of Blood 
During Nanoparticle Assisted Drug Delivery Process Through A Permeable Micro Vessel." 
Microvascular Research 92 pp 25-33, 2014 DOI:10.1016/j.mvr.2013.12.007 

[12] R Bali, N Gupta and S Mishra. “Dispersion Characteristics Of Non-Newtonian Fluid During 
Transportation Of Nanoparticles In Permeable Capillary." Applications and Applied Mathematics 
Vol 11, Issue 2, pp 632-645, 2016 

[13] J V R Reddy, D Srikanth and S K Das. “Modelling And Simulation Of Temperature And 
Concentration Dispersion In A Couple Stress Nanofluid Flow Through Stenotic Tapered Artery." 
The European Physical Journal Plus 132:365 , 2017 DOI: 10.1140/epjp/12017-11643-1 

[14] R Bali and N Gupta. “Study Of Transport Of Nanoparticles With K-L Model Through A 
Stenosed Micro Vessels." Applications and Applied Mathematics Vol 13, Issue 2, pp 1157-1170, 
2018 

[15] KM Surabhi, J V R Reddy and D Srikanth. "Impact Of Temperature And Concentration 
Dispersion On The Physiology Of Blood Nanofluid: Links To Atherosclerosis." Indian Academy of 
Sciences, Sadhana 43:210, 2018 DOI:10.1007/812406-018-0986-8 

[16] J V R Reddy and D Srikanth. “Impact Of Blood Vessel Wall Flexibility On The Temperature 
And Concentration Dispersion." Journal of Applied and Computational Mechanics. Vol 6, Issue 3, 
pp 564-581, 2020 DOI:10.22055/JACM.2019.29023.1542 

[17] S Rathore and D Srikanth. “Mathematical Study Of Transport Phenomena Of Blood 
Nanofluid In A Diseased Artery Subject To Catheterization." Indian Journal of Physics, 2021 
DOI:10.1007/s 12648-021-02166-2 

[18] SUS Choi and J A Eastman. “Enhancing Thermal Conductivity Of Fluids With Nanoparticle 
In Developments And Applications Of Non-Newtonian Flows.” International mechanical engineering 
congress and exhibition, San Francisco , U S, Conference, 1995 

[19] J Buongiorno. “Convective Transport In Nanofluids." Journal of Heat Transfer, Vol 128, Issue 
3, pp 240-250, 2006 DOI:10.1115/1.21150834 

[20] A Einstein. *Eine neue bestimmumng der molekuldimensionen." Ann Phys 19,pp 289-306, 
1906 DOI:10.1002/andp.19063240204 

[21] A Santra, S Sen and N Chakraborty. “Study Of Heat Transfer Due To Laminar Flow Of 
Copper Water Nanofluid Through Two Isothermally Heated Parallel Plates." International Journal of 
Thermal Sciences 48, pp 391-400, 2009 

[22] T Hayat, M Hussain, S A Shehzad and A Alsaedi. “Flow Of A Power Law Nanofluid Past A 
Vertical Stretching Sheet With A Convective Boundary Condition" Journal of Applied Mechanical 
and Technical Physics 57, pp 173-179, 2016 

[23] M Khan and W A Khan. “MHD Boundary Layer Flow Of A Power Law Nanofluid With New 
Flux Condition." AIP Advances 6, pp 2119-2126, 2016 

[24] S Hussain A Aziz, C M Khalique and T Aziz “Numerical Investigation Of 
Magnetohydrodynamic Slip Flow Of Power-Law Nanofluid With Temperature Dependent Viscosity 
And Thermal Conductivity Over A Permeable Surface." De Gruyter Open 15, pp 867-876, 2017 
DOL:10.1515/phys-2017-0104 

[25] Usman, A Ghaffari, I Mustafa, J Muhammad and Y Altaf. “Analysis Of Entropy Generation 
In A Power Law Nanofluid Flow Over A Stretchable Rotatory Porous Disk." Case studies in thermal 
engineering, 2021 DOI:10.1016).csite.2021.101370 


(22023, IJETMS | Impact Factor Value: 5.672 | Page 455 


International Journal of Engineering Technology and Management Sciences 
Website: ijetms.in Issue: 2 Volume No.7 March - April — 2023 
DOI:10.46647/ijetms.2023.v07102.052 ISSN: 2581-4621 


[26] El Mehdi Acchal, H Jabraoui, S Zeroual, H Loulijit, A Hasnaoui and S Ouaskit. “Modelling 
And Simulations Of Nanofluids Using Classical Molecular Dynamics: Particle Size And Temperature 
Effects On Thermal Conductivity." Advanced powder technology, 2018 
DOI:10.1016/j.apt2018.06.023 

[27] B Ankamwar. “Size And Shape Effect On Biomedical Applications Of Nanomaterials.” 
Biomedical Engineering- Technical Applications in Medicine pp 94-114, 2012 DOI:10.5772/46121 
[28] S Ganguly and S Chakraborty. “Effective Viscosity Of Nanoscale Colloidal Suspensions." 
Journal of Applied Physics, 2009 

[29] L Pasol and F Feuillibois. “Effective Viscosity Of An Inhomogeneous Dilute Suspensions 
Flowing In A Wall.” XXI International Congress of Theoretical and Applied Mechanics, Warsaw 
Poland,2004 

[30] S Shaw. “Effective Shear Augmented Dispersion Of Solutes During Nanoparticle Ssisted 
Vdrug Delivery In A Micro Vessel." The Japan society of fluid mechanics, Fluid Dynamics Research. 
.2020 DOI:10.1088/1873-7005 1/ab6617 

[31] JC Maxwell. “A Treatise On Electricity And Magnetism”. Clarendon Press, Oxford, 1873 
[32] DAG Bruggeman. “Berichnung verschiedener physikalischer Konstantein von heterogenen 
substanzen I Dielektrizitatskonstanten und leiffahigkuten der mischkorpher aus isotropen 
substanzen." Annals of Physics Vol 416, pp 636-664, 1935 DOI:10.1002/anap. 193541560705 

[33] BC Pak and Y I Cho. “Hydrodynamics And Heat Transfer Study Of Dispersed Fluids With 
Submicron Metallic Oxide Particles.” Exp heat transfer 11, pp 151-170, 1980 
DOI:10.1080/08916159808946559 

[34] Y Xuan and W Roetzel. “Conceptions For Heat Transfer Correlation Of Nanofluids.” 
International Journal of Heat Transfer 43, pp 3701-3707, 2000 DOI:10.1016/S0017-93 10(99)00369- 
5 

[35] J R Philip. “Numerical Solution Of The Diffusion Type With Diffusivity Concentration 
Dependent." Division of Plant Industry, CSIRO, Deniliquin, NSW, 1956 

[36] MG Scherberg. “The Degree Of Convergence Of A Series Of Bessel Functions." University 
of Minnesota, Minneapolis. Minn, 1932 

[37] S Ijaz and S Nadeem. “Examination Of Nanoparticles As A Drug Carrier On Blood Flow 
Through Catheterized Composite Stenosed Artery With Permeable Walls." Computer methods and 
programs in biomedicine, 133, pp 83-94, 2016 DOI:10.1016/j.cpmb.2016.05.004 

[38] J Tan, W Keller, S Sohrabi, J Yang and Y Liu. “Characterization Of Nanoparticle Dispersion 
In Red Blood Cell Suspension By The Lattice Boltzmann-Immersed Boundary Method." 
Nanomaterials, 6, 30, 2016 DOI: 10.3390/nano6020030 

[39] D Ramaya, R S Raju, A Rao and A J Chamkha. “Effects Of Velocity And Thermal Wall Slip 
On Magnetohydrodynamics Boundary Layer Viscous Flow And Heat Transfer Of A Nanofluid Over 
A Non-Linearly-Stretching Sheet: A Numerical Study.” Propulsion and power research. Volume 7, 
Issue 2, pp 182-195, 2018 DOI:10.1016/j.jppr.2018.04.003 

[40] B Ramana and G Sarojamma. “Unsteady Convective Diffusion In A Herschel-Bulkley Fluid 
In A Conduit With Interphase Mass Transfer." International journal of Mathematical modelling and 
computations. Vol 2, No.3, pp- 159-179, 2012 

[41] JW Shreffler, J E Pullan, K M Dailey, S Mallik and A E Brooks. “Overcoming Hurdles In 
Nanoparticle Clinical Translation: The Influence Of Experimental Design And Surface 
Modification." ^ International journal of molecular sciences, 20, 6056, 2019 
DOI:10.3390/ijms20236056 

[42] T Hupfeld, A Sommereyns and S Barickowski. “Analysis Of The Nanoparticle Dispersion 
And Its Effects On The Crystalline Microstructure In Carbon-Additivated PA12 Feedstock Material 
For Laser Powder Bed Fusion.” Materials, Multidisciplinary Digital Publishing Institute, 
13(15):3312, 2020 DOI: 10.3390/ma13153312 

[43] A Mukherjee and B S Mazumder. “Buoyancy Effect On Dispersion Of A Solute In A Flow 
Through A Horizontal Channel." Acta Mechanica 58, pp 137-152, 1986 DOI:10.007/BF01176596 


(22023, IJETMS | Impact Factor Value: 5.672 | Page 456 


International Journal of Engineering Technology and Management Sciences 
Website: ijetms.in Issue: 2 Volume No.7 March - April — 2023 
DOI:10.46647/ijetms.2023.v07102.052 ISSN: 2581-4621 


[44] P Bihari, M Vippola, S Schultes, M Praetner, A G Khandoga, C A Reichel, C Coester, T 
Tuomi, M Rehberg and F Krombach. “Optimized Dispersion Of Nanoparticles For Biological In 
Vitro And In Vivo Studies.” Particle and Fibre Toxicology, 5:14, 2008 DOI:10.1186/1743-8977-5- 
14 

[45]  Shwetasaibal Samanta Sahoo; Mousime Xalxo; B G Mukunda. "A Study on Tourist 
Behaviour Towards Sustainable Tourism in Karnataka". International Research Journal on 
Advanced Science Hub, 2, 5, 2020, 27-33. doi: 10.47392/irjash.2020.28 

[46] Muniyandy Elangovan; Mohamed Yousuf; Mohamed Nauman; Mohammed Nayeem. 
"Design and Development of Delivery Robot for Commercial Purpose". International Research 
Journal on Advanced Science Hub, 4, 07, 2022, 192-197. doi: 10.47392/irjash.2022.047 

[47] Manikandan N; Swaminathan G; Dinesh J; Manish Kumar S; Kishore T; Vignesh R. 
"Significant Attention in Industry and Academia for Wire Arc Additive Manufacturing (WAAM) - 
A Review". International Research Journal on Advanced Science Hub, 4, 07, 2022, 198-204. doi: 
10.47392/rjash.2022.048 

[48] Shoeb Ahmed Syed; Steve Ales; Rajesh Kumar Behera; Kamalakanta Muduli. "Challenges, 
Opportunities and Analysis of the Machining Characteristics in hybrid Aluminium Composites 
(Al6061-SiC-AI2O3 ) Produced by Stir Casting Method". International Research Journal on 
Advanced Science Hub, 4, 08, 2022, 205-216. doi: 10.47392/irjash.2022.051 

[49] Ashima Saxena; Preeti Chawla. "A Study on the Role of Demographic Variables on Online 
Payment in Delhi NCR". International Research Journal on Advanced Science Hub, 4, 08, 2022, 217- 
221. doi: 10.47392/irjash.2022.052 

[50]  Vishnupriya S; Nirsandh Ganesan; Ms. Piriyanga; Kiruthiga Devi. "Introducing Fuzzy Logic 
for Software Reliability Admeasurement". International Research Journal on Advanced Science 
Hub, 4, 09, 2022, 222-226. doi: 10.47392/irjash.2022.056 

[51] GANESAN M; Mahesh G; Baskar N. "An user friendly Scheme of Numerical Representation 
for Music Chords". International Research Journal on Advanced Science Hub, 4, 09, 2022, 227-236. 
doi: 10.47392/irjash.2022.057 

[532]  Nirsandh Ganesan; Nithya Sri Chandrasekar; Ms. Gokila; Ms. Varsha. "Decision Model 
Based Reliability Prediction Framework". International Research Journal on Advanced Science Hub, 
4, 10, 2022, 236-242. doi: 10.47392/irjash.2022.061 

[53]  Vishnupriya S; Nithya Sri Chandrasekar; Nirsandh Ganesan; Ms. Mithilaa; Ms. Jeyashree. 
"Comprehensive Analysis of Power and Handloom Market Failures and Potential Regrowth 
Options". /nternational Research Journal on Advanced Science Hub, 4, 10, 2022, 243-250. doi: 
10.47392/irjash.2022.062 

[54] Minh Duc Ly; Que Nguyen Kieu Viet. "Improvement Productivity and Quality by Using Lean 
Six Sigma: A Case Study in Mechanical Manufacturing". International Research Journal on 
Advanced Science Hub, 4, 11, 2022, 251-266. doi: 10.47392/irjash.2022.066 

[55] Ragunath A; Poonam Syal. "Net Zero Energy Buildings Initiatives - A Review". International 
Research Journal on Advanced Science Hub, 4, 11, 2022, 267-271. doi: 10.47392/irjash.2022.067 
[56] Suresh P; Justin Jayaraj K; Aravintha Prasad VC; Abishek Velavan; Mr Gokulnath. "Deep 
Learning for Covid-19 Identification: A Comparative Analysis". International Research Journal on 
Advanced Science Hub, 4, 11, 2022, 272-280. doi: 10.47392/irjash.2022.068 

[57] Chirag H B; Darshan M; Rakesh M D; Priyanka D S; Manjunath Aradya. "Prediction of 
Concrete Compressive Strength Using Artificial Neural Network". International Research Journal 
on Advanced Science Hub, 4, 11, 2022, 281-287. doi: 10.47392/irjash.2022.069 

[58] Minh Ly Duc; Que Nguyen Kieu Viet. "Analysis Affect Factors of Smart Meter A PLS-SEM 
Neural Network". International Research Journal on Advanced Science Hub, 4, 12, 2022, 288-301. 
doi: 10.47392/irjash.2022.071 

[59] Lely Novia; Muhammad Basri Wello. "Analysis of Interpersonal Skill Learning Outcomes in 
Business English Students Class". /nternational Research Journal on Advanced Science Hub, 4, 12, 
2022, 302-305. doi: 10.47392/irjash.2022.072 


(22023, IJETMS | Impact Factor Value: 5.672 | Page 457 


International Journal of Engineering Technology and Management Sciences 
Website: ijetms.in Issue: 2 Volume No.7 March - April — 2023 
DOI:10.46647/ijetms.2023.v07102.052 ISSN: 2581-4621 


[60] | Ms. Nikita; Sandeep Kumar; Prabhakar Agarwal; Manisha Bharti. "Comparison of multi-class 
motor imagery classification methods for EEG signals". International Research Journal on Advanced 
Science Hub, 4, 12, 2022, 306-311. doi: 10.47392/irjash.2022.073 

[61]  Aniket Manash; Ratan Kumar; Rakesh Kumar; Pandey S C; Saurabh Kumar. "Elastic 
properties of ferrite nanomaterials: A compilation and a review". International Research Journal on 
Advanced Science Hub, 4, 12, 2022, 312-317. doi: 10.47392/irjash.2022.074 

[62]  Prabin Kumar; Rahul Kumar; Ragul Kumar; Vivek Rai; Aniket Manash. "A Review on 
coating of steel with nanocomposite for industrial applications". International Research Journal on 
Advanced Science Hub, 4, 12, 2022, 318-323. doi: 10.47392/irjash.2022.075 

[63] Twinkle Beniwal; Vidhu K. Mathur. "Cloud Kitchens and its impact on the restaurant 
industry". International Research Journal on Advanced Science Hub, 4, 12, 2022, 324-335. doi: 
10.47392/irjash.2022.076 

[64] T. Pravin, C. Somu, R. Rajavel, M. Subramanian, P. Prince Reynold, Integrated Taguchi cum 
grey relational experimental analysis technique (GREAT) for optimization and material 
characterization of FSP surface composites on AA6061 aluminium alloys, Materials Today: 
Proceedings, Volume 33, Part 8, 2020, Pages 5156-5161, ISSN 2214-7853. 
doi.org/10.1016/j.matpr.2020.02.863. 

[65]  R. Ranjith, C. Somu, G. Tharanitharan, Venkatajalapathi.T, Naveenkumar M, Integrated 
Taguchi cum Grey Relational Experimental Analysis (GREAT) for Optimization and Machining 
Characterization of Cryogenic Cooled AA6063 Aluminium Alloys, Materials Today: Proceedings, 
Volume 18, Part 7, 2019,Pages 3597- 605, https://doi.org/10.1016/j.matpr.2019.07. 291. 

[66] R. Devi Priya, R. Sivaraj, Ajith Abraham, T. Pravin, P. Sivasankar and N. Anitha. "Multi- 
Objective Particle Swarm Optimization Based Preprocessing of Multi-Class Extremely Imbalanced 
Datasets". International Journal of Uncertainty, Fuzziness and Knowledge-Based Systems Vol. 30, 
No. 05, pp. 735-755 (2022). Dot: 10.1142/80218488522500209 

[67] M. S. N. K. Nijamudeen, G. Muthuarasu, G. Gokulkumar, A. Nagarjunan, and T. Pravin, 
"Investigation on mechanical properties of aluminium with copper and silicon carbide using powder 
metallurgy technique," Advances in Natural and Applied Sciences, vol. 11, no. 4, pp. 277—280, 2017. 
[68] Pravin T, M. Subramanian, R. Ranjith,Clarifying the phenomenon of Ultrasonic Assisted 
Electric discharge machining, “Journal of the Indian Chemical Society", Volume 99, Issue 10, 2022, 
100705, ISSN 0019-4522, Doi: 10.1016/j.jics.2022.100705 

[69] | V.S. Rajashekhar; T. Pravin; K. Thiruppathi , “Control of a snake robot with 3R joint 
mechanism", International Journal of Mechanisms and Robotic Systems (IJMRS), Vol. 4, No. 3, 
2018. Doi: 10.1504/IJMRS.2018.10017186 

[70] T. Pravin, M. Sadhasivam, and S. Raghuraman, “Optimization of process parameters of Al- 
10% Cu compacts through powder metallurgy,” Applied Mechanics and Materials, vol. 813-814, pp. 
603—607, 2010. 

[71] Rajashekhar, V., Pravin, T., Thiruppathi, K.: A review on droplet deposition manufacturing- 
a rapid prototyping technique. Int. J. Manuf. Technol Manage. 33(5), 362-383 (2019) 
https://doi.org/10.1504/IJ/[MTM.2019.103277 

[72] Rajashekhar V S, Pravin T, Thirupathi K, Raghuraman S.Modeling and Simulation of Gravity 
based Zig-zag Material Handling System for Transferring Materials in Multi Floor Industries. Indian 
Journal of Science and Technology.2015 Sep, 8(22), pp.1-6. 


(02023, IJETMS | Impact Factor Value: 5.672 | Page 458 


